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Abstract 

We consider a Markovian approximation, of weak coupling type, to an 
open system perturbation involving emission, absorption and scattering by 
reservoir quanta. The result is the general form for a quantum stochastic 
flow driven by creation, annihilation and gauge processes. A weak matrix 
limit is established for the convergence of the interaction-picture unitary 
to a unitary, adapted quantum stochastic process and of the Heisenberg 
dynamics to the corresponding quantum stochastic flow: the convergence 
strategy is similar to the quantum functional central limits introduced 
by Accardi, Frigerio and Lu'^'. The principal terms in the Dyson series 
expansions are identified and re-summed after the limit to obtain explicit 
quantum stochastic differential equations with renormalized coefficients. 
An extension of the Pule inequalities'^' allows uniform estimates for the 
Dyson series expansion for both the unitary operator and the Heisenberg 
evolution to be obtained. 

1 Introduction 



In the interaction picture, the unitary Ut arising from a time-dependent pertur- 
bation Vt, is given by 

C/t = Texp|-i^ dsVs^ (1.1) 

where T is Dyson's time-ordering operation. A principal aim of quantum field 
theory is then to obtain a normal-ordered version of Ut- When Vt involves a 
sum of monomials of canonical quantum fields, we may use Feynman rules to 
expand Ut'. we associate a vertex to each monomial, with the number of legs 
corresponding with the degree; we then construct the class ^ of Feynman di- 
agrams consisting of such vertices with certain legs contracted (internal lines) 
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and the remainder free (external lines) ; we then specify a rule for writing down 
an operator Lq (t) which, for each G S 5^, will be a normal-ordered product 
of the fields associated to the external lines of G. We then determine a de- 
velopment of the form Ut = TliGe^^Gi^)- Now, if G can be decomposed as 
two disconnected sub-diagrams Gi and G2, then Lq = N LgiLg2 where N is 
Wick's normal-ordering operation. This leads to a second presentation of Ut- 



where is the class of connected Feynman diagrams. 

If, in place of quantum fields, we considered quantum white noises, then the 
time-ordered presentation corresponds to a Stratonovich form while the normal- 
ordered presentation corresponds to an Ito form. Our aim is not to justify 
this statement, for which there is ample support [3],[1],[5], but to prove an 
asymptotic result which, effectively, is an analogue of the Wong-Zakai theorem 
for classical stochastic processes. The interaction that we shall be interested in 
is given below as (1.8), and is quadratic in the reservoir creation/annihilation 
operator fields af (A): the corresponding connected Feynman diagrams will 
have at most two legs and therefore will be linear chains. These describe a 
reservoir quanta created, subsequently multiply-scattered (i.e., at several times 
annihilated and immediately re-created) and finally reabsorbed: external lines 
may also be present. 

We shall be interested, not in the S-matrix limit t — » 00, but in the more 
subtle van Hove [6|, or weak coupling, limit where we rescale time as t/X^ 
with A a coupling strength parameter appearing in Vt and consider the limit 
A ^ with t fixed. The fields af (A) will converge, in a sense to be spelled out 
below, to quantum white noises: more correctly, integrated versions of these 
fields converge to the fundamental quantum stochastic processes of Hudson and 
Parthasarathy's theory [7]. The van Hove limit turns out to have dominant 
contribution from Feynman diagrams where there is no overlap in the time 
ranges of the individual connected subgraphs: these are the so-called type I 
terms. All other terms {type II) are suppressed. A similar feature is observed 
for the limit of the dynamical fiow of observables. 

1.1 The Classical Wong-Zakai Theorem 

Wong and Zakai [8J studied Langevin type equations driven by diffcrentiable 
noises (A) having correlation (A) (A)) = j^G {^^) which became delta- 
correlated only in the limit A — > 0. They found that the limit dynamics was 
described by a stochastic differential equation taking the same form as the pre- 
limit equations in the Stratonovich calculus. 

Let us specialize to the flow on a symplectic manifold generated by a random 
Hamiltonian 




(1.2) 




(1.3) 



a 
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where H and Fa are smooth functions on phase space and (A) are difFer- 

entiable stochastic processes converging to independent white noises. If 

is the phase trajectory starting from xq then the evolution of functions is 

Jj'*'' (/) :— f (^x^.^^^. In the hmit A ^ we obtain, in accordance with the 
Wong-Zakai resuh, the Stratonovich-Fisk equation 

dJt (•) = Jt {•, H}dt + J2 Jt {•, Fa] o dB^ (1.4) 

a 

where are independent Wiener processes and the differential is of Stratonovich 
type: here we may view the motion as that governed by the formal Hamilto- 
nian Tt = i7 + -Fq^" where ^" are white noises. A general treatment 
of these problems using the van Hove limit is well-understood [9]. These are 
the stochastic flows that preserve the Poisson bracket structure [TU]- Averag- 
ing with respect to the Wiener measure, we obtain the dynamical semigroup 
E [Jt (•)] = exp {tZ (•)}. From the Ito calculus, the generator will be the hypo- 
elliptic operator 

^{■)=Y.{{-,Fa},Fa] + {-,H} (1.5) 

Q 

which is already displayed in Hormander form. 



1.2 Quantum Markov Approximations 

It was first suggested by Spohn [TT] that the weak coupling limit should be 
properly considered as a Markovian limit underscored by a functional central 
limit. The rigorous determination of irreversible semigroup evolutions has been 
given for specific models [12],[2]- (A detailed account of the derivation of the 
master equation for a class of quantum open systems is given in Davies' book 
[13j.) The form of the generator of quantum dynamical semigroups was deduced 
[14], jT5] using the guiding principle that the semi-group be completely positive. 
Hudson and Parthasarathy [7j subsequently developed a quantum stochastic 
calculus giving an Ito theory of integration with respect to Bosonic Fock space 
processes and demonstrated how to construct dilations of the quantum dynam- 
ical semigroups mentioned above using a Fock space as auxiliary space. 

The program now is to begin with a microscopic model for a system-reservoir 
interaction and then obtain by some Markovian limit procedure, such as the 
weak coupling limit, a quantum stochastic evolution. It was first noted by von 
Waldenfels [16] that stochastic models successfully describe the weak coupling 
limit regime for the Wigner-Weisskopf atom. Later, Accardi, Frigerio and Lu 
[T] showed how to do this for an interaction of the type T^'*'-' = Eiq ® (A) H- 
£^01 ® (A) where Eiq and E^i are bounded, mutually adjoint operators on the 
system space \)s and of (A) are creation/annihilation fields having a correlation 

(a+(A)a;(A)) = lG(^^) (1.6) 
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where G (.) is integrable. In the sense of Schwartz distributions, we have 
hmA^o (or W W) = {t — s) where 7 = dt G {t) is finite. We shall 
also take an interest in the constants 

POO pO poo 

K+:= / dtG{t), K_ / dtG{t) and K := / dt \G {t)\ . (1.7) 
Jo J-00 Jo 

We shall assume that G {—t) = G (t)* so that k± = zL ia. Already in [1], 
several important steps were taken: to begin with, there is the anticipation of 
the limit algebraic structure by means of a quantum functional central limit 
theorem which captures the long time asymptotic behaviour; secondly, there is 
the identification of the principal, type I, terms in the Dyson series which survive 
the Markovian limit (they are the ones arising from only time-consecutive two- 
point contractions); finally, there is a rigorous estimate of the Dyson series 
expansion employing an argument due to Pule [2] . 

1.3 Statement of the Problem 

Our aim is to extend this result in Ij to the more general class of interactions 

t|^^ = £^11 ® a+ (A) a- (A) + ® a+ (A) + Eqi ® o^" (A) + Eoo<»l 

= E^f,C=^[4{X)]"[a;{X)f (1.8) 

( We introduce the summation convention that when the Greek indices a, (3, . . . 
are repeated then we sum each index over the values and 1 - moreover we 
understand the index a in [.]" to represent a power.) We require only the 
conditions that the system operators Eafs are bounded with ||i?ii || < 1, where 
K is the constant introduced in (1.7). 

The interaction includes a scattering term, £^11 (A) (A), and a con- 
stant term. The terms involving £01 and £10 describe the emission and ab- 
sorption of reservoir quanta and this component has been employed in models 
of laser interactions [T^ . The constant term is of little consequence as we shall 
take it to commute with the free Hamiltonian. However, the scattering term is 
highly non-trivial: we have to contend with emission, multiple scatterings and 
absorption. This means that the number of terms in the Dyson series expansion 
of 

[/(^^^ = T exp|-i dsT(^)| (1.9) 

grows rapidly (in fact, as the Bell numbers of combinatorics l21j). However, 
we are able to prove a uniform estimate of the Dyson series expansion by a 
generalization of the Pule inequalities, which we give in section 7. We are then 
able to re-sum the series to obtain an adapted, unitary process Ut of Hudson- 
Parthasarathy type (Theorem 8.1). The type of limit involved is of a weak 
character and is often referred to as convergence in matrix elements. 

We show that the Heisenberg evolution {X) = U^^'^^ {X (g) 1r) U^^^ like- 
wise converges in weak matrix elements, for fixed bounded observables X G 
B (t)s), to Jt (X) = Ut^'' {X 1r) Ut (Theorem 10.1). 
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We are able to obtain the quantum stochastic differential equations satisfied 
by Ut and by the flow Jj. In particular, these equations will involve a gauge 
differential (due to the scattering) as well as creation, annihilation and time. In 
particular, we compute the Lindblad generator for the flow. We remark that 
interactions of the type (1.8) were considered previously in the case where the 
coefficients Eap were commuting operators [TS], and Fermionic operators [12]. 
In the former strong resolvent limit was established for the common 

spectral resolution, while in the latter, the anti-commutation relations kill off 
all but type I terms. 

2 Moments and Cumulants 

Let r ([)) be the (Bose) Fock space over the one-particle Hilbcrt t). The Fock 
vacuum will be denoted by $ and the exponential vector map by e : t) i-^ F (f)). 
As usual e (0) = $. We denote the creation fields as A'^ (•), the annihilation 
fields as (•) and the differential second quantization field as dT (•), as stan- 
dard. The Weyl operator with test function / is Vl^ (/) := cxp [A+ (/) - A^ (/)] 
and we have the Weyl map W {■). 

As is well-known, the fields Q {■) = (•) + A^ (•) are Gaussian random 
fields when taken in the Fock vacuum state. More generally, we have [20] 



where H is self-adjoint on [} with spectral measure fij^ for vector state / G (}. 
This time, we are dealing with Poissonian fields. We remark that if /i^ = 
A(5i,thcn we obtain a random variable with Poisson distribution of intensity 



The coefficients S {n, m) = (—1)'^™ ^"(T) well-known combinato- 

rial factors: they are the Stirling number's of the second kind [H] and they 
count the number of ways of partitioning a set of n items into m non-empty 
subsets. 

The expansion of Poissonian field moments in terms of cumulants, or more 
generally the expansion of Green's functions in terms of their connected Green's 
functions, can best be described in the language of partitions [22j . 

A partition of the integers {1, ... ,n} is a collection of non-empty, disjoint 
subsets (called parts) whose union is {1, . . . ,n}. The set of all such partitions 
will be denoted as ^P„: there will be S {n,m) partitions of {1, . . . ,rt} having 
exactly m parts and Bn = X^m, ^ ("-j partitions of {1, . . . , n} in total. Bn are 
called the Bell numbers [^T]. 




A > 0: 



exp{A(e--l)}.^5:^^(n,m)A^ 



n m 
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Lemma (2.1) Let fi,gi, . . . , /„, Qn G t). Then 

a,/3e{0,l}" 

= n (ffi(fc)l/i(fc-l))---(ffi(3)l/i(2))(5i(2)|/i(l)> (2.1) 

^6q3„ {i(i),...,i(fc)}e^ 

where we take the various sets (parts of the partition) {i{l) , . . . ,i (k)} G A to 
be ordered so that ?' (1) < z (2) < • • • < i{k) and if the set is a singleton it is 
given the factor of unity. 

Proof. If a{i) — 0, 1, then we have the absence, respectively presence, of 
the creator (/j). Likewise /3 (i) gives the absence or presence of the i-th 
annihilator. Evidently we must have a (n) = = /? (!)• 

Essentially we have a vacuum expectation of a product of n factors (/i)]"*''^ 
[A~ (gi)]^^^^ and this ultimately when put to normal order will be a sum of terms 
each of which is a product of pair contractions {gi\fk) where i > k. For a given 
term in the sum we write i ~ k if {gi\fk) appears. An equivalence relation is 
determined by a set of contractions as follows: we always have i = i and, more 
generally, we have i = fc if there exists a sequence j {!),... j {r) such that cither 
i ^ j (1) ~ J (2) ~ . . . j (r) ~ A; or fc ~ j (1) ~ j (2) ~ ■ ■ ■ j (r) ~ i. A partition 
A in is then obtained by looking at the equivalence classes. (Singletons 
arc just the unpaired labels.) The correspondence between the terms in the 
sum and the elements of ^„ is one-to-one and the weight given to a particular 
partition A G ^„ is just the product of (5'i|/fc)'s given in (2.1). ■ 

There is a convenient diagrammatic way to understand the formula (2.1). 
We first of all associate one of four possible vertices with each component 
[A+ (/,)]"(^) \A- {g,)f'\3 = 1, ■ • • , n, they are, for (a,-, /3^.) = (1, 1) , (1, 0) , (0, 1) 
and (0, 0) respectively, 

::x r. 

Scattering Emission Absorption Neutral 
Figure 1 

We draw the n vertices in a line and proceed to join up the emission lines to 
the absorption lines (pair contractions!). A typical situation is depicted below: 
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Evidently we must again join up all creation and annihilation operators into 
pairs; we however get creation, multiple scattering and annihilation as the rule; 
otherwise we have a single neutral vertex. In the figure, we can think of a particle 
being created at vertex i (1) then scattered at i (2) ,i (3) , i (4) successively before 
being annihilated at i(5). (This component has been highlighted using thick 
lines.) Now the argument: each such component corresponds to a unique part, 
here {i (5) , * (4) , i (3) , i (2) , i (1)}, having two or more elements; singletons may 
also occur and these are just the constant term vertices. Therefore every such 
diagram corresponds uniquely to a partition of {1, . . . , n}. 

We remark that (2.1) can be considered as a special case of the expansion 

G{xi,--- ,Xn) = E^e5p„n{i(i),...,i(fe)}e.AC'(a^i(i)'-- - '^m) of an n-particle 
Green's function G in terms of the connected Green's functions C. 

Let us write ^ for the set U„^„ of finite partitions. With each partition 
A we associate a sequence of occupation numbers n = (nj)°l^ where 
Uj = 0, 1, 2, . . . counts the number of j-tuples making up A. In general, we set 

E{n):=J2:inj, N{n):=J2>ij (2.2) 

j 3 

so that if A leads to sequence n, then E (n) = n, while N (n) counts the 
number of parts making up the partition. We shall denote by the set of all 
partitions having the same occupation number sequence n. 

Given a partition A e^„ we use the convention q (j, k, r) to label the r— th 
element of the fc— th j-tuple. A simple example of a partition in *Pn is given by 
selecting in order from {1, 2, . . . , (n)} first of all ni singletons, then pairs, 
then nz triples etc. The labelling for this particular partition will be denoted 
as q {.,.,.) and explicitly we have 

q{3,k,r) = Y.^.ni + {k-l)nj+r. (2.3) 

Definition (2.2): We shall denote by 6° the collection of Pule permutations, 
that is, p e &n, E{n) = n, such that q = p° q again describes a partition in 
Specifically, 6° consists of all the permutations p for which the following 

requirements are met: 

i) the order of the individual j-tuples is preserved for each j - 

p{q{j,k,l))<p{q{j,k'A)) Vi,l<fc<fc'<n,-; (2.4) 

ii) creation always precedes annihilation in time for any contraction pair - 

p{q{j,k,l)) < p{q{j,k,2)) <■■■ < p{q{j,k,j)) Vj,l <fc<nj. 

(2.5) 
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In these notations we may rewrite the result of the lemma (2.1) as: 
Lemma (2.3) Let fi,gi, ■■■,fn,9n&^- Then 

Y: ($1 [A+ (/„)]"^"' [A- {g„)f-^ ...[A+ (A)]"^^^ [A- ig,)f'^ $) 

a,/3e{0,l}" 

E{n)=n rij j-1 

= X] "12 I[I[I[(3p{qU,k,r+l))\fp{qUM)) (2-6) 
n pee° j>2 fe=l r=l 

To better understand this, we return to our diagram conventions. Given an 
arbitrary diagram, we wish to construct the Pule permutation putting it to the 
basic form. For instance, we might have an initial segment of a diagram looking 
like the following: 



Figure 3 

There will exist a permutation a of the n vertices which will reorder the 
vertices so that we have the singletons first, then the pair contractions, then the 
triples, etc., so that we obtain a picture of the following type 

ria triples n2 pairs rii singletons 

/ixiv /ixzv ^ /zv 

Figure 4 

The permutation is again unique if we retain the induced ordering of the 
first emission times for each connected block. 

3 A Microscopic Model 

We shall consider a quantum mechanical system S (state space coupled 
to a Bose quantum field reservoir R over a one-particle space (state space 
i)n = T (t))j)). We shall take the reservoir to be in the Pock vacuum state <I>. 
The interaction between the system and the reservoir will be given by the formal 
Hamiltonian 

if(^) =Hs^lR + ls^dT (Hj,) + H^^l (3.1) 

where the operators Hs and Hj^ arc self-adjoint and bounded below on l)s and 
respectively. The interaction is taken to be 

H^!:^ = Eii^A+ {g)A- {g) + XEw®A+ {g)+XEoi(8A- {g) + X'Eoo®lR (3.2) 
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where Eap are bounded operators on i)s with En and Eqq self-adjoint and 
EiQ — Eq^. The operators (g) and (g) are the creation and annihilation 
operators with test function g G i)]^. (The parameter A is real and will later 
emerge as a rescaling parameter in which we hope to obtain a Markovian limit.) 
We shall also assume the following harmonic relations 



e " J^afS e = e Ej^fi 



-Al{g)e-"^- = 4(M). (3.3) 



where (Qt : r G M) will be the one-parameter group of unitaries on with Stone 
generator H\^. 

We transfer to the interaction picture with the help of the unitary 

U (r, A) = e+''^'^-^^'^-^"+-^^'*-^''' g-jrff(^' ^ ^^ j^ 

In the weak coupling regime, we are interested in the behaviour at long time 
scales T ~ t/}? and from our earlier specifications we see thatiJt'^^ = ^7(^/A^A) 
satisfies the interaction picture Schrodinger equation 

|c/(^) = -zT,(A)C/i"' (3.5) 

with Tf (A) as in (1.8). Here we meet the time-dependent rescaled reservoir 
fields 

a±(A):=ieT-*A'A±(0,/,.g). (3.7) 
Specifically we have 7 = dr (g, e'''^"''^-'^) g) = 2tt {g, S {H]^ - uj) g) and 



K+ = (g, —— ^- 7-r75 ) = 57 - « PV( -g ) where PV denotes 

the principle value part. 



4 Quantum Central Limit 

The limit A — *■ for the above, the two-point function becomes delta-correlated. 
However, it is vital to have a mathematical framework in which to interpret the 
limit states and observables. 
For convenience we set 

0^=exp{ir(Fj,-a;)}. (4.1) 
We assume the existence of a non-zero subspace, £, of [}]^ for which 

/oo 
\{f„d':fk)\du<^ 

whenever fj,fk G i- (In reference [1], explicit examples of dense subspaces, 
i, are given and correspond to "mass-shell" Hilbert spaces.) The question of 
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completeness can be addressed immediately: a sesquilinear form on i is defined 

by 

ifjlfk) ■■= / (/j, CA) du EE 27T {f„5 [H], - uo) h) (4.2) 



and we can quotient out the null elements for this form; the completed Hilbert 
space will again be denoted by t and (.|.) will be its inner product. The test 
vector g appearing in the interaction must belong to t so that the constant 
7 = (5 Is) is finite. 

Let W {■) be the Weyl map from t))j as before. We now fix fj E { and 
< Sj < Tj < 00 for certain indices j and introduce the rescaled operators 



At U) \ r duA^ (C/,./,) , (j) ( ^ r due-,^.fA . (4.3) 



Note that, with respect to our earlier notations (3.7), if fj — g then Af {j) = 
jj^ dua^ (A). The following result is proved as lemma 3.2 in Accardi, Frigerio 
and Lu [T]. We write 1[s,t] for the characteristic function of an interval [S, T]. 

Lemma (4.1) For the fields introduced in (4-3) 

hm [A- ij) , A+ (k)] = (/,!/,) {l[s,^T,],l[s,,m) ■ 



The right hand side is the inner product (/j (g) l[Sj,Tj]i fk ® l[s^^Tk]) on the 
Hilbert space t^L'^ {^~^)- This space is isomorphic in a natural way to the t - 
valued square-integrable functions on R+ and we denote this space as (M+ , 6) . 

The appropriate noise space for the limit A — > will in fact be the Bose Fock 
space r (i^ (R+,fi)). Indeed, we have the following fact proved as theorem 3.4 
inn. 

Theorem (4.2) Let * be the Fock vacuum for F (L^ J)) and let W (.) de- 
note the usual Weyl mapping from (R+, 5) into the unitaries on F (L^ (M+, 6)) . 
Then 

hm ($1 Wx (1) . . . Wx (fc) $) = (^1 W (A ® 1[5,,T,]) .■■W{fk® 1[5„T.]) ^) 

for arbitrary k and fj G 6 and < Sj < Tj < 00. 

5 The Dyson Series Expansion of uj:^'' 

The formal Dyson series development U^^^ — J2'^=q (~*)" (t, A) involves the 
multiple time integrals 

Dr,{t,X)= f ds„...dsiT,„(A)...T,aA). (5.1) 

JA„(t) 
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For a e ©„, we introduce the simplex 

K (*) := ■■■,si):t> s^(„) > ■ • • > 5^,(1) > 0} (5.2) 

and A„ {t) in (5.1) is the simplex corresponding to the identity permutation. 

We consider matrix elements of the type (^(j)^ W\ (1) $| u'j:^'' (j)2 (8i W\ (2) 

with G f)5 and Wa (j) as in (4.3). Substituting for the Dyson scries, wc find 
that the n-th term can be rewritten as an expectation involving the vacuum 
state $ only: 



Wx (1) $1 T«„ (A) . . . (A) </.2 ® W^A (2) 
1 $1 (A) . . . Ts, (A) (A2 (Wa (1) $1 Wa (2) $) (5.3) 



where Tg (A) is obtained from (A) by the canonical translations 

a+ (A) ^ 4 (A) + /ii (i, A) ; a," (A) ^ a," (A) + {t, A) (5.4) 

with 



That is, 



where 



/i.- (i, A) = ^ / ' (C/A^/il^VA^i/) • (5.5) 
(A) = it, A) ® [a+ (A)] " [a^ (A)] (5.6) 



Eoo{t,X) = Eoo + Eiohi{t,X)+Eoih*2{t,X) + Enhi{t,X)h*2{t,X); 

Ew{t,X) = Ew + h;{t,X)Eii; 

Eo,{t,X) = Eo,+hi{t,X)En; 

En{t,X) = En- (5.7) 

In this way we see that the n-th term in the Dyson series expansion of the 
matrix element is, up to the factor (-i)" (Wa (1) $|Wa (2) 



/ dSn---dsi ((pi\ Ea„l3^{Sn,X) . . .Eail3^{si,X) ^^2) 

X (*l [< W] K W] • • • K W] K W] (5-8) 

and our summation convention is now in place. The vacuum expectation can 
be computed using lemmas (2.1) or (2.3). The resulting terms can be split into 
two types: type I will survive the A — > limit; type II will not. They are 
distinguished as follows: 
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Type I: Terms involving contractions of time consecutive annihilator-creator pairs 
only. (That is, under the time-ordered integral in (5.8), an annihilator 
"s,+i W must be contracted with the creator af. (A).) 

Type II: All others cases. 

The terminology used here is due to Accardi, Frigerio and Lu [T]. 

We again resort to a diagrammatic convention in order to describe the Dyson 
series expansion into sums of integrals of products of two-point functions. There 
is a one-to-one correspondence between the diagrams appearing in the n-th term 
of the Dyson series and set of partitions of the n vertices. The diagram pictured 
as a typical situation in that section would contribute a weight of 

f [til, A) (ii6, A) • • • {h, A) 

"'Ai7(t) 

^Cx{tn-tii)---Cx{t2-ti) 

to the series. Let us consider a typical diagram. We shall assume that within 
the diagram there are ni singleton vertices [• • • , • • • ], '^2 contraction pairs 

[• • • ....C~^^\ -~ • • • ] , "-3 contraction triples [• • • ....C~^^^^^C~^^\ -- • • • ] , etc. That is the 
diagram has a total of n = "Y^- jrij vertices which are partitioned into m — rij 
connected subdiagrams. We see that the total number of diagrams contributing 
to the n— th level of the Dyson series will be given by the Bell number i3„. 



6 Principal Terms in the Dyson Series 

A standard technique in perturbative quantum field theory and quantum statis- 
tical mechanics is to develop a series expansion and argue on physical grounds 
that certain "principal terms" will exceed the other terms in order of magni- 
tude [231 . Often it is possible to re-sum the principal terms to obtain a useful 
representation of the dominant behaviour. Mathematically, the problem comes 
down to showing that the remaining terms are negligible in the limiting physical 
regime being considered. 

Let n be a positive integer and m e {0,...,n— 1}. Let {iPj,qj)}JLi be 
contractions pairs over indices {!,..., n} such that if P = and 
Q — {qi, ...,qm} then P and Q are both non-degenerate subsets of size m and 
we require that pj < qj for each j and that Q be ordered so that qi < ... < qm- 
We understand that {pj,qj)j^i is type I if qj = pj + 1 for each j and type 11 
otherwise. The following result is an extension of lemma 4.2 in Accardi, Frigerio 
and Lu [1] as now P DQ need not be empty. 

Lemma (6.1) Let {pj,qj)'j^i be a set of m pairs of contractions over indices 
{1, . . . , n} then 

p m 



mj.n—m 



< 



(n — my. 



(6.1) 



12 



Moreover, as A — > 0, 



« m 



0, type //. 



(6.2) 



Proof. Let q = qi and set t (q) = [5 {p) — s (q)] j)? then 



P m 

dsi...d5„ n(a;(^^.)(A)a+^^)(A) 



fds{l)...j 
Jo Jo 



s(q-2) 



ds {q - 1) 



dt 



(q) I 

Jo 



ds{q + l)... 



s{n-l) 



ds in) (^g, gj [] (^a^^^^^ (A) a+^^.) (A) 

j=2 



However, we have that s (p) — X^t (p) < s{q — 1) and so we obtain the bound 



ds {!)... 

•s(n-l) 



s(g-2) 



rfs (g - 1) / dt (q) 



s(p)-A^t(«) 



(is (? + 1) . 



/■s(,n-ij \ TT- / 

/ d5(n) (5,^r(,)5)n(«;fe)(^)<fe)W 



And so, working indiictivoly we obtain (6.1). 

Suppose now that the pairs are of type I, then p = q — 1 and so the lower 
limit of the t (g)-integral is zero. Consequently, we encounter the sequence of 
integrals 



...j^ ds{q-l)j^ dt{q)J^ ds(g+l)...(ff, 

this occurs for each g-variable and so we recognize the limit as stated in (6.2) 
for type I terms. 

If the pairs are of type II, on the other hand, then let j = min {k : pk < qk — 1}; 
setting q = qk, we encounter the sequence of integrals 



fs(q-l)/\ 



s{q-l)-\'t{q) 



fsiq-2) /. 

.../ ds{q-l) 

Jo J Is 



[s{p)-s{q-l)]/\^ 



dt{q) 



s{q-l) 



ds{q+l)...lg, 6i^(g)5()... 



but now, with respect to the variables s{l),...,s{p),...,s{q — 1) we have that, 
since s{p) ^ s{q—l), the lower limit \s [p) — s (q — 1)] /A^ of the t (g)-integral 
is ahnost always negative and so, as t {g, 9'^ g) is continuous, we have the 
dominated convergence of the whole term to zero. ■ 
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Clearly type II terms do not contribute to the n— th term in the series ex- 
pansion in the Hmit. However, we must estabUsh a uniform bound for all these 
terms when the sum over all terms is considered. We do this in the next section. 

Before proceeding let us remark that the expression (5.8) is bounded by 

Ca„/3„ . . .C„i/3i ||0i|| II02II 

X / dsr.... ds, U\ [a+ (A)] [a7„ (A)] [a+ (A)] [a" (A)] $) 

(6.3) 



where 

Cii - llSiill; 

Cw = ll^ioll + ll^iill h2\ Cm = \\Em\\ + ll^iiH h^; 

Coo = ll^ooll + ||Sio||/ii + ||£;oi||/'-2 + ||Sll||/ll/l2 (6.4) 

and h, = jr^du \{gKh)\, - J^^^du \{gKh)\. 

Recall that we require that KCu < 1 and that C = max {Cu, Cio, Cqi, Cqo} < 

oo. 

We need to do some preliminary estimation. We employ the occupation 
numbers introduced in section 2. The number of times that we will have (a, f3) = 
(1, 1) in a particular term will be J2j>2 ^ 2) nj (that is, singletons and pairs 
have none, triples have one, quadruples have two, etc.) and this equals E (n) — 
2N (n) + Til. Therefore, we shall have 

a„;3„ . . . Cc.,0^ < cf/")-2Ar(n)+„i^2W(„)-„,^ 

7 Generalized Pule Inequalities 

Putting all this together we get the bound 

Cct„/3„ • ■ • Ca^fj^ / dSn ■ ■ ■ dsi 

JA„{t) 

X (*l K K" (A)]'" . . . [< (A)]"^ [a- (A)]''^ $) 

E{n)=n 

^ ^ ^ ^S(n)-2Af(n)+ni^2Ar(n)-ni 

n peso 

X / dSn . ..dsi U U n Ca {Sp(q(j,k.r+1)) " Sp(gO,fe,r))) (7.1) 

JA„(t) -^o „_i 



where we use the estimate (6.5) and we obtain the sum over all relevant terms 
by summing over all admissible permutations of the basic q term. To estimate 
the simplicial integral we generalize an argument due to Pule (lemma 3 of [2]). 
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Let p be the induced mapping on M" obtained by permuting the Cartesian 
coordinates according to p G S^. Then the bound in (7.1) can be written as 

E{i\)=n 

^ ^S(n)-2W(n)+ni^2Ar(n)-ni 



"j J-1 



/ dSn---dsi n n n {^qU,k,r+ 



(7.2) 



j>2 fe=l r=l 



where i? = U {pA„ (t) : p £ 6^}. This is down to the fact that the image sets 
pAn (t) will be distinct for different p G 6^. Now the region, R, of integration 
is a subset of [0, i]" for which the variables Sq^j^f^ i) are ordered primarily by the 
index j and secondarily by the index k. Moreover, each of the variables 



(7.3) 



are positive, (Vj; k — 1, . . . n^; r = 1, ... ,j — 1). (These properties of R are im- 
plicit from the choice of the ordering q and of the nature of the permutations 
p € &n-) Consider the change of variables 



(Sl, . . . , Sn) Uq(j,fc,r)) 



(7.4) 



where the ordering is first by the j, second by the k, and for the it's finally by 
the r = 1, . . . , j — 1. This defines a volume-preserving map which will take R 
into A„j (t) X A„2 (t) X ■ ■ ■ X [0, oo)"^ x [0, oo)^"=* x • • • . From this we are able 
to find the upper estimate on (7.2) of the form 

^_E(n)-2Ar(n)+ni^2Ar(n)-ni (,t V IJ [t V ij 



ni 



7121 



|Ga (s)|ds 



E^_E(n)-2A'(n)+Tii^2Af(n)-ni l'^ ^ -'-J j^E(n)-N(n) 
" ni!n2!--- 



E(n)-n ^^£;(n) + BAr(n) 

rti!n2! 



""1 "n.T • • • 



where A = \n {KCn) and B = In (t V l)-|-ln (C^ V l)+ln (Cf^^ V l)-|-ln (ff-i V l) . 

The restriction to those sequences n with E (n) = n can be lifted and the 
following estimate for the entire series obtained 

^AE{n)+BN{n) °° °o ^{kA+B)nk ( ^A+B 



k=l nfc=0 



The manipulations are familiar from, for example, the calculation of the grand 
canonical partition function for the free Bose gas [24]. The requirement for 
convergence is that e"^ < 1, or equivalently, that KCn < 1. 
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8 Limit Transition Amplitudes 

We are now ready to re-sum the Dyson series. First of all, observe that the 
functions hj {t, A) defined in (5.5) will have the limits 

hj (t) := hmhj {t, A) = l[s„T,] {fj\9) • (8.1) 

Likewise, we obtain Ea/3 {t) = lim^^o -^a/S {t, A) which will be just the expres- 
sions in (5.7) with the hj {t, A) replaced by their limits. Explicitly, we have 

En [t] = En, ^oi {t) = ^ai [hi (t)]" , 

Ew (t) = Eip [hi {t)f , Eoo (t) = [hi (tT Eap [h; {t)f . (8.2) 

Secondly, only type I terms will survive the limit. This means that, for 
the n-th term in the Dyson series, the only sequences ai, 0i,a2, (^2, ■ " ,C(n,/3n 
appearing will be those for which = Q!„ = and /3; = aj+i for 1 = 1,..., n — 1. 

Thirdly, we encounter the following limit of the two point function: Gx {t — s). 
Let / and g be Schwartz functions then we will have the limit 

/ dt2 [ dtiGx{t2-ti)f{t2)g{ti)^K+[ dsf{s)g{s). 
Jo Jo Jo 

Therefore, employing lemma (2.3), we find 

lim {(f>^ ® Wx (1) $1 U (^/A^ A) (l)^ ® Wx (2) $) 

A— *0 

= (W (/i ® l[s„Til) *| W (/2 ® 1[5„T.]) *> 
~ n—1 
(-i)" / dSn--- dSi TT [K+t)+ (s;+i - Si)f' 
n •^A„(t) tl 

X ^ (<^ll^0/3„_,(s„)---^/3,/3,(s2)^/3,0(si) (8-3) 

/3e{o,i}"-' 

where we use the symbol d+ for a one-sided delta function: J d+ {t — s) f (s) ds = 

We now develop this scries. Suppose that we have /S^+i = = Pf,, that 
is, there are no contractions to the fc— th term, then we encounter the factor 
■E'OO (s) = [hi (s)]" Eai3 [/12 (.s)f where s = Sk- Otherwise, if we have contrac- 
tions on the terms associated to consecutive variables Sk+r, ■ ■ ■ , Sk+i, Sk and 
we assume that Sk+r is not paired to Sfc+r+ijUor Sk to Sk-i'. then we en- 
counter the factor Eqi {sk+r) En {sk+r-i) ■ ■ ■ En (sfe+i) Eio [sk) with the vari- 
ables Sk+r, ■ ■ ■ , Sk+i,Sk all forced equal to a common value s, say. This factor 
will then be [hi {skT E^i {Enf-^ Ei,s [h*2 {sk)f. 

Now (8.3) involves a sum over all consecutive pairings: the corresponding 
partition will have all parts consisting of consecutive labels. We can list these 
parts in increasing order, say from 1 to m if there are m of them, and let rj be 
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the size of the j-tli part . The number of contractions will be ^ /?; and this will 
be n — m = J2jLi {^j ~ !)• With these observations we see that (8.3) becomes 

ri-\ l-ry^=n 

(W (A ® 1[S„T,]) *| W (A ® 1[5„T.]) *) EE E 

[hi {SmT- [hi {Sm)f- • • • [/li {S^)r [K {Sl)f' (8-4) 

where we set 

In the following, we shall encounter the coefficients 

oo ^ 

LaH ■■= -i'^{-inY^^ E^^^p = -iEaf} - K^ai ^ _^ Eip. (8.6) 



With respect to the representation (R+, t) = ? (g) (I^"*"); we introduce 
the four fundamental operator processes (here X[o,t] is the operator on (K^) 
corresponding to multiplication by l[o,f]) 

(creation) 

(conservation) Aj^ 

(annihilation) A^^ 
(time) 

These are the basic quantum stochastic processes on the Hudson-Parthasarathy 
space r (L^ (M+, 4)) . We note that the quantum Ito table takes the concise form 

dAfdAl'^=^dAf (8.9) 
with all other pairs vanishing. 

Theorem (8.1) Suppose the system operators Eap are hounded with K < 
1. Let (pi,(p2 € and /i,/2 G 8- Then 

lim^ (^01 ® Wx (1) $1 U[^^ (f>2 (E> Wx (2) 

= ((/>i ® (/l ® l[Si,Ti]) *| C/t <^2 ® (/2 ® 1[S2,T.]) *) 

where (Ut ■ t > 0) is a unitary adapted quantum stochastic process on l)s 'Si 
r (L^ (M+,t)) satisfying the quantum stochastic differential equation 

dUt = L^fiUt ® dAf (8.10) 



= A+ (5®l[o,ti); 

= dT (\g){g\ S X[o,t]) ; 
= A- (fif0l[o,t]); 
= t. 
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with t/o = 1 and where the coefficients are given by (8.6): 

Lii = ~iEii{l + inEiiy^ , LiQ^ -i{l + iKEn)~^Eio 

Lqi = -iEoi{l +iKEiiy^, Loa = -iEoa - kEoi{1 + iKEiiy^Eia. 

Proof. The quantum stochastic differential equation (8.10) takes the form 

dUt = -{W -l)Ut(E) dA]^ + LUt ® dAl° 
7 



where 



W 1 ^ IK^Ell 

W = — — — (unitary) 

1 + iK+i^ll 

L = —i(l + «K+i?ii)^"^£'io (bounded) 

H = Ef^Q + lialn+Eai / Ew\ (self-adjoint). (8.11) 

[ I+IK+Ell J 

A fundamental result of quantum stochastic calculus 7 is that the process Ut 
defined as the solution of (8.11) with initial condition — 1, exists and is an 
adapted, unitary process. With our summation convention in place, we have 
the chaotic expansion 

^aM/3(™) . . . ^^aa)/3(l) (8 ^2) 

(t) 



^* = / ^a(m)/3(m) ' ' ' ^c(l)/3(l) ® dA'^^^)^'^'^^ ' ' ' 

and so ®W {fi® 1[Si,Ti]) *| Ut(j)2®W (/2 ^ IfSa.Ta]) can be expressed 
as 

{W (/l ® l[s„Ti]) (/2 ® 1[5„T.]) *) E ('^ll 

• • • La{2)fi{2)La{l)fj{\ 

m>0 

X [ dSra--- ds^ ([h, (S„,)]"^"^ [h^ {Sra)f"'^) ' ' ' ([hi (si)]"^'^ [/i^ {si)f' 



By inspection, this evidently agrees with (8.4). ■ 
8.1 Re-summing the Series 

Again we drop all diagrams that are type II to get the series 



o 



..Ox... 



....ax ,....+ ....cx^i.. 
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Wc sec the first appearance of scattering in the last term in the 3rd term 
of the series: such terms however eventually out-proliferate diagrams with no 
scattering. The terms have been grouped by vertex number, however, it also 
possible to group them by effective vertex number (equal to the number of parts, 
or equivalently the original simplex degree minus the number of contractions) 
to give 

O = + □ + DO + + - 

where now each box is the following sum over all effective one-vertex contribu- 
tions: 



□ 



which is analogous to the expression of the self-energy in quantum field the- 
ory:as a sum over irreducible terms. (As we have seen, one- vertex contributions 
terminate at second order when there is no scattering: as this is a form of cu- 
mulant expansion, the emission/ absorption problem is Gaussian, while allowing 
scattering means that we must have cumulant moments to all orders!) 

If the limit effective one-vertex label is t then its weight is 

-iEoo {t) + {-if kEoi {t) EiQ {t) + {-if K^Eoi (i) {t) E^ (t) + • • • 
= -iEoo {t) - kEoi {t) —^—Eio {t) 

1 -|- IK-tjii 

^ [hi{trG^0[h2{t)f 

where the geometric series can be summed since < 1. We therefore see 

that 



lim I 



(•lOel/iOlfSi.Ti])*! 



1+ / G^adA'^l^ 
Jo 



The QSDE then takes the form 

dUt = G^^dAfUt 



<S'e(/2(S)l[S2,T2])*)- 



= -{W-l)Ut® dAl^ + LUt 8) dAl^ 
7 



-L^WUt <S) dA°f - (^^iL^L + iH^ Ut O dA 
with the coefficients {W, L, H) are as before. 
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9 Dynamical Evolutions 

Let X be a bounded operator on the system state space i)s- We define its 
Heisenherg evolute to be 

) {X) := U^^'>^ [X (g) 1r] U^^\ (9.1) 

In addition, what we term the co- evolute is defined to be 

iff ^ {X) U^^^ [X ® 1r] C/f (9.2) 

We wish to study the limits of J^^^ and Kj:^^ as quantum processes taken 
relative to the Fock vacuum state ^ £ i)R for the Bose reservoir. To this end, 
we note the developments 



(9.3) 



kI^''{X) = V(-l)" / dsr,---dsiX ix)0---oX ,x^(X(E,lR). 

•^A„(t) 

(X) = V (-z)"+" / ds„ ---ds, [ dtf,--- dh 

X T(^) • • • T(^) [X (g) 1r] t|^^ • • • t[^'> , (9.4) 

where A'^/(.) 

We note that the co-evolution has the simpler form when iterated. The 
evolution itself requires a separate expansion of the unitaries. (This disparity is 
related to the proof of unitarity for quantum stochastic processes in [7], where 
the isometric property requires some work while the co-isometric property is 
established immediately.) In fact, the same inequalities as used to establish 
the convergence of uj:''^^ suffice for the co-evolution: in both cases we have a 
Picard iterated series. We remark that in [26j the co-evolution only is treated 
for emission/ absorption interactions. 

We likewise have the expansion 

01 Wx (1) $1 {X) 02 ® Wx (2) $ 



= 2, (~*)" " / ■ ■ ■ dsi / dtn ■ ■ ■ dti 

X (01 1 Ea,f3^ {si,X)...Ea„f3„ (Sn,A) X Ef, \)...E^^,, (si,A) 02) 

X ($1 [a+ (A)]"^ [a- (A)]^^ ■ • • [< (A)]"" (A)]^" 

[al (A)] K- (A)] [a+ (A)] [a" (A)] . (9.5) 

The vacuum average of the reservoir operators can be expressed as a sum of 
products of two-point functions with each summand representable as a parti- 
tion of 71 + n vertices. Our strategy is similar to before. We shall use diagrams 
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to describe the individual contributions, and attempt to obtain a uniform esti- 
mate. The Heisenberg diagrams are a more involved than last time due to the 
scattering, however, the general idea goes through again. 

Let us consider an arbitrary Heisenberg diagram. If we considered only the 
t — t contractions and ignored everything else then we would have a partition 
of the n variables, let's say with occupation numbers n = (rij). Likewise, 
if we looked at only the s — s contractions then we have a partition of the n' 
s— variables, say with occupation numbers n' = {n'j). At this stage we can then 
take the s — t contractions into account. The diagram below shows a quartet of 
s variables joined to a triple of t variables. 



s s s s t t t 

Figure 5 

Let Ijk be the number of s — t contractions joining a part of j s's to a part 
of k t's: here we use an obvious abuse of terminology, as technically they arc all 
in the same part! We also introduce the occupation numbers 1 = (Ij), 1' = (/^) 
where Ik = ^jk ^.nd = Ij^- (When no scattering was present, we only 
had the possibility that In, previously denoted as I, could be non-zero.) It is 
convenient to introduce the occupation numbers m = (m-j) and m' = ("^j) 
where nij = Uj — Ij and m'^ — rij — I'j. Here mj counts the number of parts of 
t-variables of size j having no elements contracted with an s-variable. 

The procedure adopted in the last chapter is now repeated. We consider 
equivalence classes of Heisenberg diagrams leading to the same set of sequences 
n, n', 1, 1', or equivalently m, m', 1, 1' as above. We can choose a basic Heisenberg 
diagram as the representative of each class, and there will be permutations 
p € S° and p' G ©„/ of the t and s variables respectively which will allow us 
to reorganize the basic Heisenberg diagram into any other element of the the 
class. (We omit the explicit choice of basic of Heisenberg diagram and leave its 
specification to the reader as an exercise.) 

Now for each diagram in a given class there will then be chronologically 
ordered blocks of sizes mi , m2 , • • • , m[ 

, 7712, • ■ ■ JiJ2,' ■ ■ and by the type of argument encountered before we arrive 
at the following upper bound for the sum of absolute values for all the diagrams 

E E[in+m' +l+l')-2N[m+in' +l+l')+mi+m[+h+l2 

m,m',l 

^(j2N[m+m'+l+l')-[mi+m[+li+l2) 
^^^-^^JV(m+m'+l) 



(mi!m2!---)Ki!m^!---)Gi!^2!---) 

^_E(m+m' -1-1+1') -AT (m+m'-|-l+r) Ar(l) 
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Here we add sequences of oceupation numbers componentwise, ie m + m' is 
(rrij + rn'j), etc., and we note that (1) = iV (1'). Recalling the constants A and 
B from before, and introducing S' = i In (i V 1) + In (C^ V l) + In {C^^^ V l) + 
In (^K~^ V l) + I In (7), we sum the series to get the upperbound 

which is again convergent as e"^ < 1. 

We now wish to determine the limit A — > 0. Once again, only diagrams 
having time consecutive s — s and t — t contractions, as well as non-crossing s — t 
contractions, are going to contribute to the limit. The presence of scattering 
now means that we have more diagrams, however, we can reduce this using 
the effective vertex method and, once again we can arrive at a simple recursive 
formula. This time, we have 

tO-x-Ch) = e — ) 

+ {■■■Cy x-U O ) 

Here we meet new effective vertices in the final diagram. On the right we 
have 

which for vertex time t corresponds to the operator weigth 

-iEw {t) + i-tf kEu {t) Ew {t) + {-if {t) ^11 {t) Ew {t) + ■ ■ ■ 

1 + ikEh 

= Gl0[h2{t)f. 

While on the left we have 
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which has the weight 



lEoi {t) + i'k^Eoi {t) En (<) + {K*f Eoi (t) En {t) En (*) + ••■ 

The recursion relation here is 

% O e (1) \Ul [X (g) 1] Ut(f>2 ® e (2)^ = ((/-i ® e (1) | [X ® 1] 02 ® e (2)) 
j0i®£(l)|C/4 [X® 1] (g„/3 [/Ji (ii)*]"[/i2 ih)f^dUt,q}2®£i'2) 
(01 ® £ (1) \dUl [gI^ [hi (ti)*]" [/i2 (ti)]^) [X ® 1] C/t,02 ® £ (2) 
+ / dh(^(f>,®e{l)\ul[[hi{h)XG\^XGifi[h2{h)f rgl]Ut,cj)2(S)e{2) 



A2(t) 



A2(t) 

t 



01 ® £ (1) I |l + j [(xG„/3 + G]3„X + gLxGi/s) ® 
01 ® £ (1) I |l + ^ U^Co^p (X) U.dAf^ 02 e (2)^ , 



and this is the form we want! 

To summarize, the pre-hmit flow J^"^^ : B{\)s) ^ B (hs ® hn) given by 
Jj'^'' (X) := uj;'^^^ {X (E) U^'^'' converges in the sense of weak matrix elements, 
for fixed X €B ([)s), to the limit process Jt (X) = j/ {X f ) Jt. We find that 
{Jt)t>o determines a quantum stochastic flow on f)5 F (L^ J)) and from 
the quantum stochastic calculus we obtain the quantum Langevin, or stochastic 
Heisenberg, equation 

dJt (X) ^ Jt {Cc^p {X))^dAf. 

The super-operators Cap are the required Evans-Hudson maps [27] Cap {X) = 
XGap + gI^X + G\„XGif3 and these can be written in the standard form 



Cn {X) ^ - {W^XW - X) ■ 

Cw{X) = W^[X,L\- Coi{X) = -[X,L^]W- 



Coo (X) 



2 [L\x]L + lL^[X,L]-t[X,H]. 



In particular, £oo is a generator of Lindblad type [TS]. We shall give a more 
detailed treatment of the convergence in the next section. 
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10 The Convergence of the Heisenberg Evolu- 
tion 



We now wish to determine the hmit A ^ of (9.5). We have an integration 
over a double simplex rc^gion and the main features emerge from examining 
the vacuum expectation of the product of creation and annihilation operators. 
Evidently, the vacuum expectation can be decomposed as a sum over products 
of two point functions and it is here that lemma 6.1 becomes important. What 
must happen for a term to survive the limit? If we have any contractions 
between vertices labelled by the t's then the term will vanish if the times are 
not consecutive. The same is true for contractions between vertices labelled by 
the s's. From our estimate in the previous section, we can ignore the terms that 
do not comply with this. 

As a result, contractions between the s's, say, will come in time-consecutive 
blocks: for instance, we will typically have m blocks of sizes ri, r2, • • • , (these 
are integers 1,2,3,..., and Yl^i fj = n). With a similar situation for the t's, we 
obtain the expansion 



1 ® Wx (1) $1 Jr {X) 02 ® Wx (2) $ 



X! X! H dSn---dSi dtn-- 



■dtl 



m '3 m -3 



j=i fe=i 



j=i fe=i 



x($| 



a+ (A) 



«„(i) (A) 



a^,^^ (A) 



a^(^) (A) 



+ negligible terms 

where we relabel the times as 

JJ) . 



Pi 



(A) 



a (™) (A) 



<i) (A) 

h 



a (1) (A) 



(10.1) 



and introduce the block product of system operators 

• • • 4f ; a) := E^.,, (4'^ a) (4^^ a) (sgL,; a) E,, (4f 
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Wc now examine the limit of (10.1). The estimate on the series expansion of 
the Heisenberg evolute given in the previous section shows that we can ignore 
the so-called negligible terms in (10.1). The limit is rather difficult to see at 
this stage. However, what we can do is to recast the expression that we claim 
will be the limit, 

(.^1 ® (/i ® l[s„T,]) M/| Jt {X) ct>^®W {h ® 1[S„T.]) *) , (10.2) 

with Jt {X) = Ul {X (g) 1) Ut, in a more explicit form. 

Recall the chaotic expansion of the process Ut given in (8.12), the expression 
(10.2) then becomes 

E/ / E E (--)'"^"'" 

X {W (A ® l[5„T,i) *| dA^^-P- ■ ■ ■ dA^^f^dAtf^ ■ . . dA'i^''' W {h ® l[s„T.]) *> . 

Now the expectation between the states W {fj ® l[g^ j^j) 'I' can be converted 
into an expectation between the Fock vacuum state \E' if we make the following 
replacements 

d^ii ^ dA^"" + hldA^° + hidA^^ + hihldA^° 

^ (10.3) 
where hj {t) = 1[Sj,Tj] {fj\9) as in (8.1). This leads to the development 

E / / EE 

X {<t>i I Cj, (si) ■ . . (5„.) X &^2L (i-) ■ • ■ E^ill ^2) 

X (*| rfA^-'^- • • • dA'^l^^dA'^^"'^ ■ ■ ■ dA'il''^ . (10.4) 

where the operators £'''1' {t) are given by 



Ki{t)(^Eu{t)y ^Ei0{t), r>2 



Again we note that the operators Ea/s {t) have been introduced in (8.2). 

It remains to be shown that the limit of (10.1) will be (10.4). We observe 
that 
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lim (^(j)^ ® Wa (1) $1 ^ {X) </.2 ® Wx (2) 
= ^ ^ ^ X! / • • • dsi / dijji • • • dti 

n,n m.mri.-,r„;i,...,i^-^A„(t) JA^(t) 

X lim (<i>| [a+ (A)] - K- (A)]''^---[<(A)]"'" [a-JA)]^" 

A — *U 

[< (A)] [«r. (A)] "'^ . . . [< (A)] [ar, (A)] $) . (10.5) 
We now require the fact that 

hm / ds„, ■ ■ ■ dsidt^, ■■■dh{<i>\ [a+ (A)] [a" (A)] ■ ■ ■ [a+^ (A)] [a'^ (A)] 

[al^ (A)] (A)] ''"^ . . . [a+ (A)] ''^ [a,- (A)] $) / (s,„, ■ ■ ■ , s,,t^, ■ ■ ■ .h) 

= [ {^\dA':-^- ■ ■ ■ dA':^^^dA';^^'''^ ■ ■ ■ dA';;"' ^) / (s™,--- ,si,t™,--- ,ti) 
Jr 

for / continuous and R a bounded region in m+m dimensions which is the union 
of simphces of the type (5.2). This is readily seen, of course, by expanding the 
^-expectation as a sum of products of two-point functions and reassembhng the 
Hmit in terms of the ^'-expectations of the processes A"^ . This is evident from 
theorems 4.2 and 6.1 quoted earher and from the quantum Ito calculus [?'. 

We therefore see that the limit form as given in (10.5) agrees with the stated 
limit. 

Theorem (10.1) Suppose that Eap are bounded with K \\Eii\\ < 1, as before. 
Let 01, 02 G and /i, /2 G t. Then, for X e B (t)s), 

lim^ <j0i O Wx (1) $1 j/^' {X) 02 «) Wx (2) $^ 

= (01 ® (/l ® l[Si.Ti]) *| Jt (X) 02 ® (/2 ® 1[S2,T.]) • 

To summarize, the pre-limit flow jj:^^ : B (l)s) ^ B ([)s ® \)r) given by 
Jj^-* (X) := u[^'^^ {X ® In) Uf'^'^ converges in the sense of weak matrix elements, 
for fixed X B ([}s), to the hmit process Jt {X) = u} {X ®\)Uf We find that 
{Jt)t>Q determines a quantum stochastic flow on tjs^r (L^ (M+, t)) and from the 
quantum stochastic calculus |7j we obtain the quantum Langevin, or stochastic 
Heisenberg, equation 

dJt {X) = Jt {L^p {X)) ® dAf (10.6) 

where 

The analogous result will hold for the co-evolution. Though, as mentioned 
before, there is a more immediate proof using the original estimates. 
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11 Conclusions 



We began with a discussion of time-ordered versus normal ordered presentations 
of unitary operators relating to scattering dynamics. It is suggestive to write 
the limit unitary Ut as either 



Ut 



or Ut 



T exp 



N exp 



I I ds E, 



t 

dsLap 





■a/3 



[air [a/ 



(11.1a) 
(11.1b) 



Here af are just symbols and we mean nothing more than that (11.16) is the 
solution of (8.10) while (11.1a) reminds us that we have the limit generated by 

a perturbation = Eap ® a\ (A) [at (A)]''. (Formally, of course, we might 

consider af as a limiting form of the fields af (A).) 

Remarkably, these identifications (11.1a, 6) can be viewed as presentations 
of (1.1) and (1.2) if we supplement the operators with the following white 
noise CCR 

[at, at] = {t- s) + [t - s) (11.2) 

where 5± are future/past delta functions: we would have the right hand side 
7i5 (t — s) if it was not for the fact that we invariably meet with simplicial 
integrals. The stochastic Schrodinger equation (8.10) can be written as 



LapUt [atf dt 



If we understand that 



(11.3) 

^ X^p (t) [atf dt 



dUt = 

which is in normal ordered form 

means Xap {t)'S)dA'^^ then we recover the Hudson-Parthasarathy calculus. The 
product of two quantum stochastic integrals will have to be put into normal 
order, using (11.2), but this will be equivalent to the usual quantum Ito rule 
with Ito table (8.9). 

Alternatively, we could consider the equation dUt — —iTtUt dt with T( = 

Eap ® 0.1 [cit] : this is what is suggested by (11.16). (The Hamiltonian Tj 
plays an analogous role to the one encountered earlier for classical stochastic 
Hamiltonian flows leading to (1.4).) However, the expression TtUt contains 
terms like atUt which are out of normal order and so cannot be directly inter- 
preted in the quantum Ito calculus. Nevertheless, the following purely formal 
manipulations can be used [5] 



[at.Ut] = 



at, 1 - i / TsUsds 



Eip [asf^+ {t-s)Usds 



-iK+EiiatUt 



-EioUt, 



leading to 



atUt = 



1 



1 + iK+Eii 



{Utat -iK+EioUt}. 



(11.4) 
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(Similar manipulations have been performed separately for emission-absorption 
and for scattering interactions in [28j.) 

By making the replacement (11.4), wherever it occurs, we obtain a proper 
normal ordered form and this turns out to be precisely (11.3). In the classical 
problem for the limit of the flow under the Hamiltonian (1.3), the canonical 
structure is never lost - though we have to look to the Stratonovich calculus to 
see it. We similarly have that the canonical structure is retained in the quantum 
problem - and we even have a formal Hamiltonian Tt- provided that we look at 
things in the appropriate way. 

Acknowledgement 1 The author is greatful to Ramon van Handel for many 
stimulating discussions about the original paper that lead to several improve- 
ments and a revision of the Heisenberg flow convergence proof. 
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